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Phase transition in the Sznajd model with independence. 
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Abstract. - We propose a model of opinion dynamics which describes two major types of social 
influence - conformity and independence. Conformity in our model is described by the so called 
outflow dynamics (known as Sznajd model). According to sociologists' suggestions, we introduce 
also a second type of social influence, known in social psychology as independence. Various social 
experiments have shown that the level of conformity depends on the society. We introduce this 
level as a parameter of the model and show that there is a continuous phase transition between 
conformity and independence. 



Introduction. — Opinion dynamics is one of the most 
studied subjects in the field of sociophysics. Among a 
number of models that have been proposed (for a recent 
review see pQ), simple models based on Ising spin variables 
are particularly interesting. In all these models, from the 
voter model [2], majority rule [US] to the Sznajd model 
[5], opinions are described by discrete variables S = ±1. 
The ferromagnetic state is an attractor for all three models 
in one and two dimensions, as well as in the case of com- 
plete graphs [T] . Obviously, in real social systems complete 
unanimity is never reached. Moreover, in real systems the 
public opinion does not reach any fixed point and per- 
manently changes. To make models of opinion dynamics 
more realistic, Galam has proposed two modifications: 

• Contrarian behavior [5] - with a certain probability 
an agent adopts the choice opposite to the prevailing 
choice of the others, whatever this choice is. 

• Inflexibles [7] - inflexible agents keep their opinion 
always unchanged. 

As shown by Galam, for a low concentration of contrarians 
a new mixed phase is stabilized, with a coexistence of both 
opinions, i.e. minority persists. Moreover, there is a phase 
transition into a new disordered phase with no dominat- 
ing opinion. In the case of the Sznajd model, contrarian 
behavior has been studied for the first time by de la Lama 
et al. [5] and the same results have been obtained. As 
suggested by Galam [6] , these results may be put in paral- 
lel with "hung elections" in America (2000) and Germany 



(2002). 

In the field of social psychology, contrarian behavior, 
introduced by Galam in [5], is nothing more than anti- 
conformity - a particular type of nonconformity [9] . There 
are two widely recognized types of nonconformity: anti- 
conformity and independence. From a social point of view, 
it is very important to distinguish between independence 
and anti-conformity [10] . The term 'independence' imply- 
ing the failure of attempted group influence. Independent 
individuals evaluate situations independently of the group 
norm. On the contrary, anti-conformists are similar to 
conformers in the sense that both take cognizance of the 
group norm - conformers agree with the norm, anticon- 
formers disagree. As noticed in |10j : This behaviour is a 
bit of a paradox, because in order to be vigilant about not 
doing what is expected, one must always be aware of what 
is expected. In contrast, truly independent people are obliv- 
ious to what is expected. Numerous studies have shown 
that the level of conformity /anti-conformity depends on 
the society (culture, age, etc.) [TOlITT] . 

In this paper we examine the influence of independence 
on the Sznajd model in one and two dimensions, as well 
as on a complete graph. We show that contrarians are not 
needed in order to obtain phase transitions as in [6, 8 and 
similar results can be observed in the presence of inde- 
pendent behavior. Below a certain critical independence 
value p c , the minority opinion coexists with the major- 
ity and above p c there is no majority in the society, the 
system is in the so called stalemate state. 
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To generalize our model and make it more realistic from 
a social point of view, we introduce also a flexibility factor 
/, which describes what is the probability of an opinion 
change in the case of independent behavior. We show, 
both analytically (for complete graphs) and using Monte 
Carlo simulations, that the critical threshold of indepen- 
dence p c decays with the flexibility factor /. This means 
that in an inflexible (conservative) society the critical in- 
dependence factor is high. This result implicates that in 
conservative societies, even if the level of independence is 
high, there is always a majority in the system - in the case 
of democratic voting, one of the two options wins. On the 
contrary, in highly flexible societies the stalemate situation 
is highly probable. Our result is one of the possible ex- 
planations why in modern societies (in which "tradition" 
is not as valuable as in old days), status quo situations 
happen more and more often. 

Model. — We consider a set of N individuals, which 
are described by the binary variables: S = 1 (f) or S = — 1 
(i). At each elementary time step, a group of people is 
chosen randomly and it influences its surrounding indi- 
viduals. In the original model [5] only one type of social 
influence (conformity) was considered: 

1. On a complete graph, two individuals are chosen at 
random and they influence a third randomly chosen 
individual |12j . 

2. In one dimension (ID), a pair of neighboring individu- 
als SiSi-\-\ is chosen and it influences two neighboring 
sites Si— i, Si+2- In this paper, to be consistent with 
the case of a complete graph, we will use the modified 
version in which only one of the two (left Si_i or right 
Si-fa), chosen randomly, will be changed. This kind 
of modification has been introduced for the first time 
by Slanina [15] . 

3. Several possibilities of generalization to the square 
lattice were proposed by Stauffer et al [14] . Here we 
use probably the most popular rule - a 2 x 2 panel of 
four neighbors is chosen randomly and influences its 
surroundings. In this paper we use a modified version 
to be consistent with the rules above - only one from 
the 8 neighbors of the panel is randomly chosen to be 
changed. 

In real social systems, conformity is one of the most 
recognized types of social response. There are many fac- 
tors that affect the likelihood of conformity, among them 
culture is one of the most important. For example, the 
level of conformity (and simultaneously anti-conformity) 
is much higher in Japan than in America [10] , The well 
known American slogan "Do your own thing" reflects a 
tendency to independent behavior. Therefore, in this pa- 
per we introduce a second type of social response, known 
as independence. With probability p, an individual Sk 
chosen to be changed will not follow the group, but act in- 
dependently - with probability / it will flip, i.e. Sk — > —Sk 



and with probability 1 — / stay unchanged, i.e Sk — > Sk- 
The parameter / is called flexibility, since it describes how 
often an individual will change its opinion in case of inde- 
pendent behavior. With probability 1 — p the individual 
will follow the usual Sznajd conformity rules, described 
above. 

Model on a Complete Graph. We consider a set 
of N Ising spins Si = ±1, i = 1, . . . ,N on a complete 
graph. In each elementary time step t two spins S, and Sj 
are chosen randomly. They will influence a third randomly 
chosen spin Sk in the following way: 

• conformity (original Sznajd rule), with probability 1 — 
p: S k {t + dt) = Si(t) if Si(t) = Sj(t), otherwise S k (t + 
dt) = S fe (t), 

• independence, with probability p: Sk(t + dt) = 
—Sk(t), with probability / or Sk(t + dt) = Sk{t), with 
probability 1 — /. 

Time t — > t + 1 after N elementary time steps, i.e. dt 
IV. 

In the case of a complete graph, the state of the system 
is completely described by the magnetization (or public 
opinion from a social point of view) defined as: 



1 N 

(*) = ]v£*(*) 



(i) 



Let us denote by Nf the number of spins 'up' and by 
iVj, the number of spins 'down': 



Nf-N^ 



N 

mN. 



(2) 



From the equations above we can easily derive the formula 
for the probability of choosing randomly a spin S = + 1 : 



P + (t) = ^= 1 + WW 
+ W N 2 



(3) 



The evolution of the system is described by the following 
equation: 

Pi-P+ =p/P_-p/P + +(l-p)P?P_-(l-p)p2p+, (4) 



where we use the notation P' + = P. 



(t + dt) and P + = 
P+(t) and the probability of choosing a 'down' spin is 
P_ = I — P + . We look for fixed points of the above trans- 
formation: 



P+-P+ = (l-2P+)[(l-p)i^-(l-p)P + 
The equation above has the following solutions: 

Pi = 1/2 for p€ [0,1] 

li2 l-p± v /(A) - 
Pi. = — — — — f° r p < 
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Fig. 1: Stationary value of magnetization for two values of the 
flexibility factor, / = 0.5 (solid line) and / = 0.25 (dotted line), 
from analytical calculations and from Monte Carlo simulations 
for the system size TV = 10 3 . Averaging was done over 10 3 
samples. 



where A = (1— f>)(l— p— 4p/). Fixed points for the magne- 
tization can be easily calculated from relation ([3]) and are 
presented in Figure [TJ There is a continuous phase tran- 
sition at p c = 1/(1 + 4/) - for p < p c minority coexists 
with majority and for p > p c there is a stalemate (status- 
quo) situation. To confirm our analytical results, we have 
provided Monte Carlo simulations on a complete graph, 
for several lattice sizes. Initially the system has been fer- 
romagnetically ordered (m(0) = 1) and then evolved ac- 
cording to the algorithm described in this section. These 
results are presented in Figure [T] and agree with the ana- 
lytical prediction. 

Scaling. — Although the model considered in this pa- 
per has 2 parameters, independence p and flexibility /, 
it can be shown that in fact results depend on the ratio 
pf /(l —p + pf) (see Figure[2]). To show this scaling let us 
first sum up the rules of the model as: 

Rule I the usual Sznajd rule, 

Rule II a site is fliped (Si — > —Si), 

Rule III nothing happens. 

At each iteration, we choose a group of sites that will 
attempt to convince one of their neighbours, chosen at 
random, then one of the 3 rules above is chosen at random 
to be followed by this chosen neighbour, with the rules 
I, II and III being followed with probability 1 — p, pf 
and p(l — /) respectively. Firstly, we note that rule III 
doesn't change the state of the system (the opinions of 
each of the sites) and that the parameters p and / are 
not important in a given iteration, once the rule to be 
followed is chosen. Suppose now that we keep records of 
the state of the system after the iterations where either 
rule I or rule II was followed. These records would form 



a sequence of states, whose statistical properties depend 
only on the initial conditions and on the ratio between the 
probabilities of following rule I and II: 



Pf 



1 



P 



(7) 



So if two models have parameters (p, /) and (p',f) such 
that r = r', then they would both generate statistically 
similar sequences, given the same initial conditions. In 
order to compare both models, we need to find the mean 
time for a change r (measured in iterations) as a func- 
tion of p and /. We follow either rule I or rule II with 
probability p = 1 — p + pf, so: 



t = p + 2p(l - p) + 3p(l - p) 2 + ... 
= p((l + (l-p) + (l-p) 2 + ...) 
+ ((l-p) + (l-p) 2 + (l-p) 3 + 
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Therefore: 



r(pj) 



1 



l-p+pf 



0) 



To compare the models with parameters (p, f) and (p' , f) 
we must look at time instants t and t' , such that 



f 



t(pJ) r(p',f) 



(10) 



Finally, we are interested in the behavior of the magne- 
tization m(t,p, f) of the model with parameters (p, f) at 
the time instant t. From 



t(l-p + P f) = t'(l-p'+p'f) 

pf/a-p) = p'/'/(i-p') 



follows: 



m{t,p,f)=m(t',p',f). 
If we pick /' = 1, then 

t' = t(l-p + pf) 



P 



So defining M(t,p) 



m{t,p,f) = Mlt{l-p + pf) 



: Pf /(l-p+pf) 

m(t,p, 1) it follows that 
Pf 



1-P + Pf 



and for each value of t, plotting 

m(t/(i-p+ p f),p,f) x Pf /(i-p+pf) 
should collapse the curves, as 

Pf 



t 

1-P + Pf 



,pj 



M t 



I-P + Pf, 



(11) 
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Fig. 2: Rescaled results from Fig. [T] 



We have used the scaling given by Eq. (|T3[) and indeed the 
data collapses in the case of a complete graph (see FigJ^J) . 
In the next sections we will see that the scaling derived 
here is valid also in the case of one- and two-dimensional 
lattices. 

In figure [2] we can see a phase transition between the 
situation where a majority exists and the stalemate situa- 
tion. This transition happens for a value of independence 
p ci that depends on the flexibility value /. We can apply 
this scaling to find that this dependence must be of the 
form 



1 



Pc 



1 + a/' 



(17) 



where, if P c is the critical independence for / = 1, then 
a is given by 



a 



1 

T- 1 - 



(18) 



This means that if there is a phase transition, then the 
critical value of independence decreases with increasing 
flexibility. We will give some interpretations for this fur- 
ther on, when the model in a two-dimensional lattice is 
analyzed. 

Model on a one-dimensional lattice. We con- 
sider a chain of length ./V with periodic boundary condi- 
tions. Each site i = 1, . . . , N of the chain is occupied by 
an Ising spin Si = ±1. At each time step we choose ran- 
domly a spin Si and side s (s — 1 for right, s — — 1 for 
left.) The updated state is: 

• conformity (original Sznajd rule), with probability 1 — 
p: Si(t + dt) = S l+S (t) if S l+S (t) = S l+ 2s(t), otherwise 
Si{t + dt) =Si(t), 

• independence, with probability p: Si(t+dt) — —Si(t), 
with probability / or Si (t + dt) — Si (t) , with proba- 
bilityl-/. 




Fig. 3: Dependence between the magnetization m and the in- 
dependence factor p for several values of the flexibility fac- 
tor /, for a one dimensional system of size TV = 10 4 (upper 
panel). The system evolved ('termalized') from an initially or- 
dered state, m(0) = 1. 'Termalization' time was r = 500MCS 
and averaging was done over 10 samples. The scaling found 
in the case of a complete graph still works (bottom panel). 



Time t — > t + 1 after N elementary time steps, i.e. 
dt = 1/N. We have chosen as an initial condition, fer- 
romagnetic order (m(0) = 1). We made Monte Carlo sim- 
ulations for several lattice sizes N — 4 x 10 2 , 9 x 10 2 , 64 x 
10 2 , 10 4 ,4 x 10 4 , but here we present results for only one 
selected lattice size, N — 10 4 . We have measured the 
magnetization of the system after various 'termalization' 
times t £ [10, 10 4 ]. The averaging has been done over 10 3 
samples. First of all, we notice that the scaling found in 
the previous section and given by Eq. (fT5|) is still valid for 
one dimensional systems - see Fig. [3] 

The results presented in Fig. |3]suggest the existence of a 
phase transition, which is unexpected for one-dimensional 
systems with short-range interactions. However, one 
should remember that the initial state in our simulations 
was ordered (m(0) = 1). Starting from this ordered state, 
the system evolves toward a stationary state that depends 
on the ratio pf/(l —p + pf)- After a short 'termalization' 
time, like in Fig. [31 the system might be still ordered but 
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Fig. 4: Dependence between the magnetization m and the in- 
dependence factor p, for a flexibility factor / = 0.5, several 
'termalization' times, r, and an one dimensional system of size 
TV = 10 4 . The system evolved ('termalized') from an initially 
ordered state, m = 1. Averaging was done over 10 3 samples. 
With an increasing r, the critical value p c is decreasing, sug- 
gesting the lack of the phase transition in one dimension. 



finally it will reach its real steady state, which is expected 
to be disordered. To check the validity of our expectations 
let us now present the dependence between the magnetiza- 
tion m and the independence factor p, for several 'termal- 
ization' times r and a given value of flexibility / = 1/2. It 
has been seen that with an increasing r the threshold value 
p* , below which the system is ordered, decreases, suggest- 
ing the lack of the phase transition in one dimension. For 
the infinite system, N — > oo and r — > oo, order is present 
in the system only for pf/(l — p + pf) = 0. This is an 
expected result, because only short range interactions are 
present in the model and the independence pf/ (1 —p+pf) 
plays the role of a temperature. 

Model on a square lattice. In this section we 
consider a square lattice L x L with periodic boundary 
conditions. Each lattice site is occupied by an individual, 
characterized by a binary opinion Si = +1 (in favor) or 
Si = — 1 (against). In each elementary time step t, a 
2x2 box of four neighboring spins is chosen randomly 
and influences one of the 8 neighboring sites of the box, 
denoted as Sf. 

• • 

• t t • 

• ft- 



In this paper, we use a modified version of the two- 
dimensional model introduced to study duopoly markets 
in [TS], therefore we update the state in the following way: 

• conformity, with probability 1 — p: If all four spins in 
the box have the same value, they will convince one 



of the eight nearest neighbors Si, changing its orien- 
tation in the direction of the spins in the box. If one 
of the spins in the box has the opposite orientation to 
the other three spins, then the neighbor changes its 
orientation to the orientation of the majority, with 
probability 3/4. In the case when there is no major- 
ity, i.e., two spins in the panel are up and two are 
down, nothing changes, i.e. Si(t + dt) = Si(t). 

• independence, with probability p: Si(t+dt) — —Si(t), 
with probability / or Si(t + dt) = Si(t), with proba- 
bilityl-/. 

We measure again the magnetization as a function of the 
independence factor p for several values of the flexibility 
factor / (see Fig. [SJ . In the case of the square lattice there 
is no doubt that there is a well defined continuous phase 
transition (see Fig. [5]) at p = p c . As in the other cases, 
scaling (given by eq. (H3J)) is valid and the critical value 
p c of independence increases with decreasing flexibility /. 

This means that in the case of high / (non-conservative 
societies) consensus is possible only in the case of low in- 
dependence (high conformity). For example (see Fig. [5]), 
for a relatively conservative society with / = 1/16 the 
critical independence is p c ~ 0.8. This means that for a 
level of independence up to 4/5 there is a majority in the 
society. If the society is less 'conservative', e.g. / = 1/4 
then the critical value of independence is p c « 0.5. This 
means that if the level of independence were, for example, 
p = 0.7 there would be no majority in the system (status- 
quo), while for the same value of p in the case of / = 1/16 
consensus would be possible. This result could be one of 
the possible explanations for why the status-quo situation 
is more and more common in modern societies. Another 
explanation, connected with the public debate, has been 
proposed several years ago by Galam [6] . 

Summary. — In this paper we have introduced a 
modified version of the original Sznajd model, in which 
two types of social influence were considered - indepen- 
dence (with probability p) and conformity (with proba- 
bility 1 — p). Conformity in our model has been mod- 
eled analogously to the Sznajd model, i.e. in the case of 
conformal behavior, individuals have followed the group 
norm. Of course, this type of social influence could be 
modeled also by the voter or majority models. In the case 
of independent behavior individuals take actions (change 
opinion) independently on the group norm. Of course, 
even in the case of independent behavior an individual 
can change opinion, but it does not depend on the social 
norm. Therefore, we have introduced a flexibility factor 
/, which denotes the probability of opinion changes in the 
case of independent behavior. For / = 0, individuals never 
change opinion in the case of independent behavior and for 
/ = 1, the opinion is always changed. From this point of 
view, varying / we can model the level of conservatism in 
the society - it decreases with increasing /. 

We have studied the model in three cases: complete 
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Fig. 5: Dependence between the magnetization m and the in- 
dependence factor p for several values of the flexibility factor /, 
in the case of a square lattice, 101 x 101. The system evolved 
('termalized') from an initially ordered state, m(0) = 1. 'Ter- 
malization' time was r = 500MCS and averaging was done 
over 10 3 samples. The critical value of independence p c decays 
with flexibility / (upper panel). The scaling found in the case 
of a complete graph still works (bottom panel). 



graphs, one-dimensional systems and two-dimensional 
square lattices. We found that, in the case of the complete 
graph and of the two-dimensional system, there is phase 
transition for a critical value of independence p c — jzf^j > 
where a is a constant that depends on the lattice. In the 
case of a complete graph it can be calculated analytically 
as a = 4. Below the critical value of independence, p < p c , 
the majority coexists with the minority, i.e. the public 
opinion is m ^ 0. Therefore, in the case of a democratic 
voting, one of the two options wins. For high indepen- 
dence, p > Pc, there is a stalemate situation in the society, 
i.e. m = 0. In our opinion the most interesting result is 
the dependence between the critical value of independence 
and flexibility, p c — j^j ~ the critical value of indepen- 
dence decreasing with increasing flexibility /. 

This means that in the case of high / (non-conservative 
societies) consensus is possible only in the case of low in- 
dependence (high conformity). On the other hand, in con- 




Fig. 6: Dependence between the magnetization m and the in- 
dependence factor p, for a flexibility factor / = 0.5, several 
'termalization' times, r, and a square lattice, 101 x 101. The 
system evolved ('termalized') from an initially ordered state, 
m(0) = 1. Averaging was done over 10 3 samples. 



servative societies, even in the case of high independence, 
consensus is possible. In modern societies the value of 
tradition and hence the level of conservatism seems to be 
decreasing. This, according to our model, could be one of 
the possible explanations for why the status-quo situation 
observed by Galam [5] is more and more encountered now 
days. 
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